Astronomy & Astrophysics manuscript no. lendl'AA201219585 


©ESO 2012 


August 2, 2012 





WASP-42 b and WASP-49 b: two new transiting sub-Jupiters. * ** 

M. Lendl 1 , D. R. Anderson 2 , A. Collier-Cameron 3 , A. P. Doyle 2 , M. Gillon 4 , C. Hellier 2 , E. Jehin 4 , T. A. Lister 5 , 
P. F. L. Maxted 2 , F. Pepe 1 , D. Pollacco 6 , D. Queloz 1 , B. Smalley 2 , D. Segransan 1 , A. M. S. Smith 2 , A. H. M. J. Triaud 1 , 

S. Udry 1 , R. G. West 7 , and P. J. Wheatley 8 

1 Observatoire de Geneve, Universite de Geneve, Chemin des maillettes 51, 1290 Sauverny, Switzerland e-mail: 
monika . lendl@unige . ch 

2 Astrophysics Group, Keele University, Staffordshire, ST5 5BG, United Kingdom 

3 School of Physics and Astronomy, University of St. Andrews, North Haugh, Fife, KY16 9SS, United Kingdom 

4 Universite de Liege, Allee du 6 aout 17, Sart Tilman, Liege 1, Belgium 

5 Las Cumbres Observatory, 6740 Cortona Drive Suite 102, Goleta, CA 93 1 17, USA 

6 Astrophysics Research Centre, School of Mathematics & Physics, Queen's University, University Road, Belfast BT7 INN, United 
Kingdom 

7 Department of Physics and Astronomy, University of Leicester, Leicester, LEI 7RH, United Kingdom 

8 Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom 



ABSTRACT 

We report the discovery of two new transiting planets from the WASP survey. WASP-42 b is a 0.500 ± 0.035 Mj planet orbiting a Kl 
star at a separation of 0.0548 ±0.0017 AU with a period of 4.98 16872 ± 7.3 x 10~ 6 days. The radius of WASP-42 b is 1.080 ±0.057 Rj 
while its equilibrium temperature is T eq = 995 ± 34 K. We detect some evidence for a small but non-zero eccentricity of e = 
0.060 ± 0.013. WASP-49 b is a 0.378 ± 0.027 Mj planet around an old G6 star. It has a period of 2.7817387 ± 5.6 x 10~ 6 days and a 
separation of 0.0379 ± 0.001 1 AU. This planet is slightly bloated, having a radius of 1.1 15 ± 0.047 R j and an equilibrium temperature 
of T eq = 1369 ± 39 K. Both planets have been followed up photometrically, and in total we have obtained 5 full and one partial transit 
light curves of WASP-42 and 4 full and one partial light curves of WASP-49 using the Euler-Swiss, TRAPPIST and Faulkes South 
telescopes. 

Key words, binaries: eclipsing - planetary systems - stars: individual: WASP-42 - stars: individual: WASP-49 - techniques: spec- 
troscopic - techniques: photometric 



1. Introduction 

Since the discovery of the first extrasolar planet around a Solar- 
type star by |Mayor & Queloz (1995]) over 700 exoplanets have 
been found. In recent years an increasingly large number of plan- 
ets have been discovered by transit surveys, i.e. surveys that 
search for planets that pass in front of their host stars producing 
characteristic photometric signals. At the time of writing, the 
two most prominent ground-based transit surveys are HATnet 
dBakos et al.|2004| > and WASP ( |PoUacco et al.|2006) . 

From ground-based transit surveys, a population of hot 
close-in giant planets has been explored. These Hot Jupiters 
are located at small separations, having periods of typically 2 
to 5 days. While the formation of gas giants is accepted to 
take place in cooler regions of the protoplanetary disc at sev- 
eral AU from the star, the process by which the planets have 
migrated to the proximity of the star has been a matter of de- 
bate. Migration can occur from interactions with the disc itself 



* Based on photometric observations made with WASP-South, 
EulerCam on Euler-Swiss telescope, the Belgian TRAPPIST telescope, 
the Faulkes South Telescope and spectroscopic observations obtained 
with CORALIE on the Euler-Swiss telescope and HARPS on the ESO 
3.6 m telescope (Prog. ID: 087.C-0649) 

** The photometric time series and radial velocity data in this work 
are only available in electronic form at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi- 
|bin/qcat?J/A+A7| 



(Goldreich & Tremaine 1980] Lin et al. T996J), or due to dynam- 
ical interactions between massive bodies in the system, namely 
scattering between pl anets (|Rasio & Ford|1996 1 and Kozai mi- 
gration (Kozai 1962; Eggleton & Kiseleva-Eggleton 2001 Wu 
& Murray 2003), which put the planet on a highly eccentric or- 
bit that then becomes tidally circularized. Recent measurements 
of the sky-projected spin-orbit angles have revealed several mis- 
aligned systems ( |Hebrard et a l. 2008; Triau d et al.|2010 Winn 
et al.|2010) , a fact which could be caused by planets undergoing 
Kozai migratio n ( Fabry c ky & Tre maine 200~7| |Wu et al.j| |2007; 
Triaud et al.|2010) . 



Many of these planets have been shown to posses radii which 
are larger than expected from models of irradiated planets (e.g. 
Fortney et al. 2007 ). Low densities appear to be common for 



anets (e.g. WASP-39 b,|Faedi et al 



201 1 J> and up 



lartman et al. 



Saturn-mass p 

to planets of the mass of Jupiter (e.g. HAT-P-32 b 
201 1 ). The most striking examples are WASP-17 b having a den 
sity of p = 0.06 pj ( |Anderson et al.|2010|[20TT]l and Kepler- 12 b 
with a density of p = 0.09 pj ( [Fortney et al.||2011 1. It appears 
that there is a mechanism depositing energy into the planet and 
thus inflating it or slowing the contraction of the planet since 
its formation. A variety of mechanisms have been proposed to 
account for this effect: the deposition of kinetic energy stem- 
ming from strong winds driven by large day/night temperature 
contrasts (Showman & Guillot 2002); enhanced opacities due to 
higher planetary metallicity causing a slow-down in contraction 
( |Burrows et al.|2007] ); heating by tidal forces from the circular- 
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Fig. 1. Phase-folded WASP photometry of WASP-42 (top) and 
WASP-49 (bottom). The unbinned data are depicted in gray; the 
black points represent the same data set binned to 2 minutes. 



ization, synchronization and re-alignment of the planetary orbit 
( Bodenheimer et al.|2001 1; reduced heat transport efficiency by 
layered convection inside the planet ( Chabrier & Baraffe|2007 1, 
and Ohmic heating from currents induced through winds in the 
planetary atmosphere (Baty gm"& Stevenson|2010 1. It has been 
found that the degree to which the planet is bloated is correlated 

Enoch 



with the incident stellar flux (Demory & Seager 



20TT1 



et al.|20lT 2012 Laughlin et al.|201 1) , favouring models which 
incorporate the absorption of stellar radiation as a cause of the 
planetary bloating. 

In this paper, we announce the discovery of two additional 
transiting planets from the WASP-South survey. WASP-42 b is a 
0.5 Mj planet in a 5 day orbit around a Kl star, and WASP-49 b, 
a bloated 0.4 Mj planet, is orbiting a metal poor G6 star every 
2.8 days. In Section [2] we present the discovery and follow-up 
observations of WASP-42 and WASP-49, leading to their iden- 
tification as harbouring transiting extrasolar planets. In Section 
[3] we describe the analysis of our data before putting the two 
planets into context in Section |4] 

2. Observations 

2.1. WASP -Photometry 

The objects WASP-42 (2MASS 12515557^204249) and 
WASP-49 (2MASS 06042146-1657550) were observed using 
the WASP survey telescopes. The WASP survey is operated 
from two sites, one in each hemisphere: the Observatorio del 
Roque de los Muchachos in the Canary Islands in the North, 
and the Sutherland Station of the South African Astronomical 
Observatory (SAAO) in the South. Each site is equipped with 
eight commercial 1 1 cm, f=200 mm Canon lenses on a single 
mount. Details of the WASP survey, its photometric reduction 
and candidate selection process can be found in |Pollacco et al. 
(2006 1 and Collier Cameron et al. (2007 ). Both targets presented 
in this work were observed exclusively from the southern WASP 
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Fig. 2. Upper panel: Radial Velocity observations of WASP-42 
together with the circular (dashed line) and eccentric (solid line) 
models. Black dots denote CORALIE data, while blue trian- 
gles denote HARPS data. Middle panel: residuals of the circular 
model. Lower panel: residuals of the eccentric model. 



site. In total 25880 data points were obtained for WASP-42 be- 
tween May 2006 and April 2008, while 18461 data points were 
obtained for WASP-49 between October 2006 and March 2010. 
The data obtained for both targets are shown in Figure [T] In both 
cases a periodically occurring transit-like signal was detected in 
the data using the transit-search algorithms of Collier Cameron 
et al. (2006 1, triggering a closer inspection of the targets and 
their selection for further study. 



2.2. Spectroscopic Observations 

WASP-42 and WASP-49 were observed with the CORALIE 
spectrograph mounted on the 1.2 m Euler-Swiss telescope at the 
La Silla site (Chile). For WASP-42 and WASP-49, 27 and 25 
data points were obtained from April 2010 to March 2011, and 
from August 2009 to October 2011, respectively. For all spec- 
troscopic observations we obtained radial velocities using the 



weighted Cross-Correlation technique (Ba ranne et al.|[T996] l as 
implemented in the CORALIE and HARPS reduction pipelines. 
The radial velocities folded on the transit ephemeris are shown 
in Figures[2]and[3] In order to ensure that the radial-velocity vari- 
ations are not caused by star spots we checked the CCF bisector 
spans (shown in Figure|4| according to the method described by 
Queloz et al. (2001 ). After the planetary nature of the signal had 



been confirmed, WASP-42 was also observed with the HARPS 
spectrograph (Mayor et al. 12003) which is located at the ESO 



3.6 m telescope at La Silla observatory. 35 data points were ob- 
tained with HARPS, 25 of which were observed nearly consec- 
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(a) Bisector measurements of WASP-42 against the measured ra- 
dial velocity. Black dots denote data points from CORALIE while 
blue triangles denote points obtained with HARPS. 

Fig. 4. Bisector spans measured for WASP-42 (left) and WASP- 
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Fig. 3. Radial Velocity observations of WASP-49 together with 
the respective model and residuals. All data have been obtained 
with CORALIE. 



utively on 4 April 2011 (UT) in order to measure the Rossiter- 
McLaughlin effect ( jRossiter|1924| |McLaughlin||1924| l and thus 
determine the projected spin-orbit angle of the WASP-42 sys- 
tem. However, due to an inaccurate ephemeris at the time of 
planning of these observations, the continuous observations be- 
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(b) Bisector measurements of WASP-49 against the measured ra- 
dial velocity. All data have been obtained with CORALIE. 



(right). 

gan only during the egress of the planet (see Figure |6| and thus 
the Rossiter-McLaughlin effect was not measured. 

2.3. Follow-up Photometry 

Two transits each of both WASP-42 b and WASP-49 b were 
observed using EulerCam between January and March 2011. 
EulerCam is an e2v 4k x 4k back-illuminated deep-depletion sil- 
icon CCD detector which was installed at the Cassegrain focus 
of the 1.2 m Euler-Swiss telescope in September 2010. The field 
of view of EulerCam is 15.68 x 15.73 arcmin, producing a res- 
olution of 0.23 arcsec per pixel. The instrument is operated at a 
temperature of -1 15 °C, cool enough to make dark current negli- 
gible (measured at less than 2.6 e~ in 30 minutes). EulerCam can 
be read out using either one or four ports, giving readout times 
of 6.5 s and 25 s. For single-port readout the readout noise is 
~ 5 e~/pixel. For four-port readout the readout noise is the same 
except that the lower-left port shows a slightly elevated noise of 
~ 8 e~/pixel. As flat-field uncertainties are often a limiting factor 
in high-precision photometry, EulerCam uses a feedback scheme 
for the telescope guiding to keep the stars on the same locations 
on the detector during the observations. This system, dubbed 
"Absolute Tracking" is based on a combination of the SCAMP 
([Berlin 2006 ) and Sextractor (Bert in & Arnouts|1996| l software 
packages. After recording an image, the positions of the stars 
on the image are extracted and then matched with a catalogue. 
From this match, the offset of the telescope from the nominal 
position is calculated using a PID algorithm and the telescope 
pointing is adjusted between exposures in order to compensate 
for drifts. Due to technical problems the 4-port readout mode 
was not available during the first transit of WASP-49, leading to 
a slightly degraded time sampling. All EulerCam observations 
were done with a defocus of 0. 1 mm in order to improve the duty 
cycle and spread the light over more pixels and thereby improve 
the sampling of the PSF 



M. Lendl et al.: WASP-42 b and WASP-49 b: two new transiting sub-Jupiters. 





0.70 



_J I I I I , I , I , I I I I , I I L_ 

-0.05 0.00 0.05 0.10 
Time since mid-tronsit [d] 



0.15 



-0.05 0.00 0.05 

Time since mid-tronsit [d] 



Fig. 5. Follow-up photometry of WASP-42 (left) and WASP-49 f right) together with the respective model and residuals. Each light 
curve has been divided by the respective photometric baseline model as described in Section 3.2 The light curves of WASP-42 
4re (from top to bottom) EulerCam on March 20 and 25 2011, and TRAPPIST on March 5, 10, 30 and April 4 2011 (UT). The 
light curves of WASP-49 are (from top to bottom) EulerCam on January 19 and March 24 201 1, TRAPPIST on January 19 and 24 
October 201 1, and FTS on 21 March 201 1 (UT). 
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Four transits of WASP-42 b and two transits of WASP-49 b 
were observed with the automated Belgian TRAPPIST tele- 
scope, also located at La Silla. For details of TRAPPIST see 
|Gillon et aL] ( |2011| > and |Jehin et al.| ( |2011) . Again, the telescope 
was defocused giving a FWHM of 3.2 arcsec on the images. 

Finally, one partial transit of WASP-49 b was observed with 
the Faulkes Telescope South (FTS) which is based at Siding 
Spring, Australia. Details of all observations are summarized in 
Table [T] and the light curves are depicted in Figure [5] 

All follow-up light curves were obtained from bias- and 
flat-field-corrected images using relative Aperture Photometry 
where several apertures were tested and reference stars were 
chosen with great care. In the case of TRAPPIST, IRAFQ was 
used in the reduction process. For EulerCam, the standard reduc- 
tion procedure is as follows. After overscan, bias and flat-field 
correction, the photometry is extracted from the images for sev- 
eral circular apertures with radii ranging from 12 to 50 pixels 
placed on the target and all other bright stars in the field. The 
optimal combination of reference stars is found by using the tar- 
get itself to measure the quality of the references. We start the 
process with the reference star with which we obtain the lowest 
RMS transit light curve and iteratively add those other references 
which yield the best improvement at any given step. Typically, 
the final reference sources is made of between three and ten sta- 
ble stars of similar brightness and color to the target and which 
are not affected by any short-term variations (caused, e.g., by 
proximity to bad pixels). 

During the observations of WASP-49 with EulerCam we no- 
ticed a faint source, 9.2 arcsec North of the target. It is outside 
the smaller apertures for which we extracted photometry in the 
analysis of the EulerCam data. As the transit is present in light 
curves created using these apertures (while the data are nois- 
ier than for larger apertures), we are certain the transit is on the 
main target. However, in order to evaluate the contamination for 
larger apertures, we re-observed WASP-49 on 15 January 2012 
with EulerCam without defocusing the telescope. The observa- 
tions were performed outside of transit, a total of 8 images were 
obtained over 10 minutes using an r' Gunn filter. We find the 
brightness ratio between the target and the contaminant to be 
1 : 0.00343 + (5 x 10~ 5 ), and, as its impact on the observed tran- 
sit depth is much smaller (4 X 10~ 5 ) than the l-<x-error quoted in 
Table [3] we neglect it in the analysis presented below. 



3. Determination of System Parameters 

3.1. Stellar Parameters 
3.1.1. Rotation period 

We have analyzed the WASP lightcurves of WASP-42 and 
WASP-49 to determine whether they show periodic modula- 
tion due to the combination of magnetic activity and the ro- 
tation of the star. We used the sine-wave fitting method de- 
scribed in Maxted et al. ( 201 l| l to calculate periodograms for 
both stars over 4096 uniformly spaced frequencies from to 2.5 
cycles/day. The false alarm probability (FAP) for the strongest 
peak in these periodograms was calculated using the boot-strap 
Monte Carlo method also described in Maxte d et al.| ( |201 l| l. 
Variability due to star spots is not expected to be coherent on 
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Fig. 6. Photometric and radial velocity observations obtained 
on April 4 2011 (UT). Upper panel: TRAPPIST transit light 
curve binned in 2 minute intervals. Lower panel: the the HARPS 
time series observation aiming at measuring the Rossiter- 
McLaughlin effect. Depicted are the residual values from a 
purely Keplerian model. The dashed red line indicates the time 
of fourth contact. 



long timescales as a consequence of the finite lifetime of star- 
spots and differential rotation in the photosphere, and so for both 
stars we analyzed the data from each observing season indepen- 
dently. 

We did not find any significant periodic signals (FAP< 
0.05) in our data apart from frequencies near 1 cycle/day and 
its harmonics. We examined the distribution of amplitudes for 
the most significant frequency in each Monte Carlo trial and 
used these results to estimate a 95 % upper confidence limit of 
1 milli-magnitude for the amplitude of any periodic signal in the 
lightcurves for both WASP-42 and WASP-49. 



3.1.2. Spectroscopic Analysis 

For each star a total of 1 1 (WASP-42) and 23 (WASP-49) in- 
dividual CORALIE spectra were co-added to produce single 
spectra with typical S/N values of 60:1 (WASP-42) and 100:1 
(WASP-49). The standard pipeline reduction products were used 
in the analysis. 

The analyses were performed using the methods given in 
IGillon et al.| ( |2009l ). The H a line was used to determine the ef- 
fective temperature (7^), while the Na i D and Mg i b lines 
were used as surface gravity (logg) diagnostics. The parameters 
obtained from the analysis are listed in Table [2] The elemental 
abundances were determined from equivalent width measure- 
ments of several clean and unblended lines. Values for micro- 
turbulence (£ t ) were determined from Fe i using the method of 
Magain (1984). The quoted error estimates include those given 
by the uncertainties in T s g, log g and as well as the scatter due 
to measurement and atomic data uncertainties. 

The projected stellar rotation velocities (vsinz„) were de- 
termined by fitting the profiles of several unblended Fe i lines. 
Values for macroturbulence (v mac ) of 1.4 + 0.3 kms -1 (WASP- 
42) and 2.9 + 0.3 kins"' (WASP-49) were assumed, based on 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



the tabulation by Gray (2008). For both cases, an instrumental 
FWHM of 0. 1 1 ± 0.01 A was determined from the telluric lines 
around 6300A. Best fitting values of vsin;* = 2.7 + 0.4 kms~' 
(WASP-42) and vsim* = 0.9 + 0.3 kms" 1 (WASP-49) were 
obtained. 

There is no significant detection of lithium in the spectra 
of either star, with equivalent width upper limits of 4 mA, cor- 
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target 


date (UT) 


telescope/instrument 


Filter 


T r \HJD T Rn - 


2450000] 


Pre/ 


RMS [relative flux, per 2 min] 


WA <iP AO 


os rn 901 1 


TP APPTCT 
1 Knrrlo 1 




JOZj.OJOZl ± 


OOOzlA 


1.1 J 


1 9 v 1 O -3 
1 .Z A 1U 


WA CP AO 


i o o^ on 1 1 


TP APPTCT 
1 KArrlo 1 


ItZ 


C/C^n AAOQS -i- 

joju.O'+uyo ± 


o ooo^ 


1 .4 J 


1 A v 1 n -3 


WA CP AO 


?n rn 901 1 

ZU-Uj-ZUI 1 


EulerCam 


r Gunn 


SAzLO ^099 1 -i 


- OOAO 


1 .tO 


O Q v 1 O -3 

\j.y x iu 


WA CP AO 


9S 0^ 901 1 


EulerCam 


r Gunn 


S^AS SS^OO -+- 
OCHO.joOUU ± 


O 000^7 


1 -UD 


1 Q v 1 O -3 

1 .7 X 1U 


WA CP AO 


an O^ 901 1 


TP APPTCT 
1 KArrlo 1 




JDJU. DO 1 1 1 ± 


o ooo^ 


1 1 8 
1.10 


1 a v 1 n -3 

1 .^f X 1 u 


WASP-42 


04-04-2011 


TRAPPIST 


I+ z ' 


5655 54921 + 


00056 


1.15 


1.2 x 10~ 3 


WASP-49 


19-01-2011 


EulerCam 


r' Gunn 


5580.59412 ± 


0.00041 


1.26 


0.5 x 10" 3 


WASP-49 


19-01-2011 


TRAPPIST 


I+z 1 


5580.59523 ± 


0.00074 


1.00 


1.5 x 10" 3 


WASP-49 


21-02-2011 


FTS 


z' Gunn 


5613.976 ± 


0.011 


1.00 


0.7 x 10" 3 


WASP-49 


24-03-2011 


EulerCam 


r' Gunn 


5644.57521 ± 


0.00035 


1.00 


0.7 x 10" 3 


WASP-49 


24-10-2011 


TRAPPIST 


I+z' 


5858.76832 ± 


0.00048 


1.08 


1.4 X 10" 3 



Table 1. Summary of follow-up photometry. Target, date, telescope and filter are given for each observation together with the mid- 
transit time, red noise amplitude a nd the RMS of the binned (2 minutes) residuals. 
(1) as defined inlWinn et al.|(|2010 1 



responding to an abundance upper limit of log A(Li) < 0.5 + 
0.2 (WASP-42) and 8mA, corresponding to an abundance upper 
limit of log A(Li) < 0.7 ± 0. 1 (WASP-49). 



Parameter 


WASP-42 


WASP-49 


RA 


12 h 51 m 55.62 s 


06 h 04 m 21.47 s 


DEC 


-42°04'25.2" 


-16°57'55.1" 


V mag 


12.57 


11.36 




5200 ± 150 K 


5600 ± 150 K 


log 8 


4.5 ±0.1 


4.5 ±0.1 


ft 


0.8 ± 0.2 km s" 1 


0.9 ±0.2 km s~' 


v sin i. 


2.7 ± 0.4 km s"' 


0.9 ±0.3 km s~' 


[Fe/H] 


0.05 ±0.13 


-0.23 ± 0.07 


[Na/H] 


0.23 ±0.10 


-0.10 ±0.08 


[Mg/H] 


0.21 ± 0.10 


-0.09 ± 0.06 


[Si/H] 


0.18 ±0.06 


-0.04 ± 0.05 


[Ca/H] 


0.11 ± 0.12 


-0.06 ±0.11 


[Sc/H] 


0.13 ±0.09 


0.06 ± 0.05 


[Ti/H] 


0.27 ± 0.15 


0.05 ± 0.06 


[V/H] 


0.59 ±0.13 


0.02 ±0.16 


[Cr/H] 


0.12 ± 0.12 


-0.16 ± 0.03 


[Co/H] 


0.30 ± 0.06 


-0.03 ± 0.09 


mm 


0.14 ±0.09 


-0.14 ±0.06 


log A(Li) 


< 0.5 ± 0.2 


< 0.7 ±0.1 


Mass 


0.89 ± 0.08 M Q 


0.94 ± 0.07 M Q 


Radius 


0.87 ±0.11% 


0.90 ±0.11% 


Sp. Type 


Kl 


G6 


Distance 


160 ± 40 pc 


170 ± 20 pc 



Table 2. Stellar parameters of WASP-42 and WASP-49 from 
Spectroscopic Analysis. 

Note: Mass and Radius estimates using the Torres et al. ( 2010| > 
calibration. Spectral Type estimated from T e g using the table in 
|Grayl ( [2008] l. 



3.2. Combined Analysis 

The combined analysis of discovery and follow-up photome- 
try and spectroscopic data was performed in several steps with 
the Markov Chain Monte Carlo (MCMC) code described in 
|Gillon et aL] pOTUl [20T2] > The code makes use o f the tran- 
sit light curve models by |Mandel & Agol| ((2002 ) and uses 
a Keplerian model to fit the radial- velocity measurements. 
For radial-velocity measurements obtained during transit, the 
prescription of the Rossiter-McLaughlin effect provided by 



Gimenez ( 2006[ ) is used. Various models for the photometric 
baseline (e.g. 0th to 4th order polynomials with respect to time 
and external variables such as pixel shifts) can be included in 
the fit of the transit light curves via minimization of the model 
coefficients at each MCMC step. The basic jump parameters are 
transit depth dF , impact parameter b, transit duration d, epoch of 
mid-transit Tq, period P and Ki — K Vl - e 2 P x ^ (where K and 
e denote the radial velocity semi-amplitude and eccentricity, re- 
spectively). For these parameters we assume a uniform prior dis- 
tribution. In order to take into account the stellar limb-darkening, 
we adopt a quadratic law interpolating the coefficients tabulated 
by |Claret & Bloemen ( 2011| l. However, we chose their combi- 
nation c\ = 2u\ + U2 and C2 = u\ — 2u2 as jump parameters as 



introduced by |Holman et al. ( 2006| l. The physical basis of these 
parameters leads us to assume a normal prior distribution with a 
width equal to the errors of these parameters. Where eccentric- 
ity is not set to e — 0, it is included in the analysis via including 
the jump parameters yfe cos w and yfe sin u> (where cj denotes 
the argument of periastron). We use the calibration technique 
devised by [Enoch et al.| ( |2010[ ) which uses a fit to a set of well 
studied main-sequence eclipsing binaries in order to infer the 
stellar mass and radius from the mean stellar density (measured 
directly from the transit light curves, [Seager & Mallen-Ornelas| 
(2003), temperature and metallicity. This technique is based on 
a similar relation by Torres et al. ( 2010| > where the stellar surface 
gravity is used in the place of the stellar density. A minimum of 
two chains is run in order to check their convergence using the 
Gelman & Rubin test ( |Gelman & Rubin|1992| l. 

As a first step, the photometric data obtained by the WASP- 
South cameras were analyzed with the aim of finding a mean 
mid-transit time and period estimate from these long-term data. 
As the light curves from WASP contain many large outliers (i.e. 
points deviating by more than 10 times the transit depth), we 
discarded points deviating by more than 5cr before running two 
MCMC chains of 10 4 points each. Here the only parameters left 
to vary were the time of mid-transit and the period, while all 
other parameters were kept fixed at their approximate values. 
The period derived here was used as a starting point in the sub- 
sequent analysis of the high-precision follow-up data, while the 
mid-transit time was included as an extra constraint on the pe- 
riod. 

Then, for each follow-up lightcurve we performed MCMC 
analyses including only the light curve itself and the CORALIE 
radial velocity measurements. We tested various models for the 
photometric baseline, using chains of 10 5 points. A second-order 
polynomial in time was assumed as a minimal baseline model to 
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(a) Modified Hertzsprung-Russell diagram showing the location 
of WASP-42 together with isochrones for ages of (from bottom 
to top) 0.5, 1.0, 2.0, 3.0, 5.0, 10.0 and 13.0 Gyr. The red dashed 
lines indicate evolutionary paths for (from left to right) 1.0, 0.9 and 
0.8 Mq. Isochrones and evolutionary paths have been interpolated 
from |Marigo et al.| ( [2008] > using z = 0.021. 
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6500 6000 5500 5000 
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(b) Modified Hertzsprung-Russell diagram showing the location 
of WASP-49 together with isochrones for ages of (from bottom to 
top) 0.5, 1.0, 2.0, 3.0, 5.0, 10.0 and 13.0 Gyr. The red dashed lines 
indicate evolutionary paths for (from left to right) 1.1, 1.0, 0.9 and 
0.8 Mq. Isochrones and evolutionary paths have been interpolated 
from |Marigo et al.|p008) using z = 0.01 1 . 



Fig. 7. The locations of WASP-42 and WASP-49 in the modified Hertzsprung-Russell diagram. 



take care of stellar variability, airmass and other time-dependent 
effects. More complicated models (i.e. polynomials up to the 
fourth degree with respect to time, position shifts, sky back- 
ground, and FWHM) were tested with the Bayesian Information 
Criterion ( Schwarz 1978 ). Only for the EulerCam light curve ob- 
tained on 19 January 201 1 (WASP-49) was a more complicated 
model (a linear fit to position plus the quadratic fit to time) found 
to be justified. 

Finally, we performed a global analysis in order to find a 
single solution using all available follow-up data together with 
the mid-transit time obtained from the discovery photometry. 
The dependencies and orders in the photometric baseline models 
were set to those found before. For each planet we ran three sets 
of two chains of 10 5 points. From the first set, we determined 



discussed in Section|4]while the planetary and stellar parameters 
of both objects are presented in Table[3] 



4. Discussion 

We announce the discovery of two new transiting extrasolar 
planets in the southern hemisphere from the WASP-South sur- 
vey. 

4.1. WASP-42 b 



WASP-42 b is a 0.5 Mj planet orbiting the Kl star 
2MASS 12515557^1204249 every 5 days. The host-star metal- 
licity is near the solar value. 

Based on the v sin z, and stellar radius of WASP-42, we de- 
duce a maximal rotation period of P = 16.1+3.2 d. Using 
et al.| ( |2008) . In the further analysis, the photometric errors were the |Barnes| ( |2007| l relation, this gives an upper age limit of 



correction factors for our photometric errors based on the calcu- 
lation of the red noise ( |Pont et aT7||2Q06| > as described in |Winn| 



multiplied by these factors. For the RV data, we determined "jit- 
ter" factors which serve to scale the measurement errors to the 
standard deviation of the best-fit model. These factors are added 
quadratically to the RV errors and compensate both instrumen- 
tal and astrophysical effects (such as stellar activity) that are not 
included in the initial error calculation. The values for the RV 
jitter are m s" 1 (WASP-42, CORALIE), 2.4 m s _1 (WASP-42, 
HARPS) and 10 m s' 1 (WASP-49, CORALIE). 

The second and third sets of chains were run with the adapted 
errors, first setting the eccentricity to e — and then leaving 
it free. For WASP-49, no significant deviation from a circular 
orbit is detected (e = O.OIS^q qjo), and all parameters agree well 
within one sigma for circular and eccentric case. The eccentric 
solution for WASP-42 yields e = 0.060+° ^ and will be further 



0.84 + n« Gy from gyrochronology. We derived a second es- 



-0.35 

timate based on the R', 



HK activity indices (|Noyes et al.||1984 
from the HARPS spectra following the procedures ot|Lovis et 



84) 

[aL] 



( [20T1) . We find \og( R' HK ) = -4.9 + 0.7, indicative of a rather 



quiet star. Using the Mamajek & Hillenbrand (2008! relation 
we derive a rotation period of P = 40.1 + 4.7 d and an age 
of 6.1 + 1.2 Gy. We also interpolated the isochrones of Marigo 



et al. ( 2008| l using z = 0.021 and compared them to the location 



of WASP-42 in the modified ((p*/Po) vs - T es) Hertzsprung- 



Russell diagram (Figure 7(a) I. Here we used the temperature 



determined from the spectroscopic analysis together with the 
stellar density determined from the global analysis described in 
Section 3.2 It should be noted that in the analysis of the spectra 
a degeneracy exists between macroturbulence v mac and stellar ro- 
tation v sin i m . As the value for v mae is not directly measured but 
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WASP-42 


WASP-49 




e = 




e = 


Jump parameters 


AF = {R p IR,f 

V = a* cos(; p ) [R,] 

T u [d] 

T m - 2450000 [HJD] 
P[d] 

K 2 [kms-'a 113 ] 

Cl/ 

C2f 
Cl,I+z' 
t'2,I+z' 
C'l.z' 

C2,z' 


0.01650 ± 0.00039 
0410 +o.o45 

U.tiu _ .055 

0.12042 + S« 

-0.00095 

5650.56723 ± 0.00024 

A QS 1 0000068 

4-98168 //_o 0000073 
109.8 ± 2.9 
1.271 ±0.056 
0.122 ±0.063 
0.908 ± 50.00 
-0.234 ± 0.033 


0.01650 ±0.00037 
0.4183$ 
0.12043 ± 0.00093 
5650.56720 ± 0.00023 
4.9816872 ± 0.0000073 
110.4 ±2.9 
1.270 ±0.055 
0.120 ±0.065 
0.907 ± -0.049 
-0.234 ± 0.032 


0.01376 ±0.00038 

0.745 ± 0.014 

0.08832 ± 0.00080 

5580.59436 ± 0.00029 

2.7817387 ±0.0000056 

79.9 ±3.3 
i 174+0.058 

1 ' 1 ^-0.061 

-0.130 ±0.048 
0.835 ± 0.043 
-0.241 ±0.046 
0.786 ± 0.042 
-0.257 ± 0.027 


Deduced parameters 


K [km s" 1 ] 

R P [RA 

M„ [MA 

e 

Oi [deg] 
a [AU] 
a/R, 
i P [deg] 

b„ 

T occ - 2450000 [HJD]" 

P P IpA 
T eq [K]* 

M» [Mq] 

R, [Sq] 

P. \PQ] 
UiJ 

U2,i> 

«l,I+z' 

«2,I+z' 

«l,z' 

«2,z' 


64.3 ± 1.7 
1.063 ±0.051 
0.497 ± 0.035 


0.0547 ± 0.0017 
13.84 ± 0.34 

0410 +«:645 

5658.03977 ± 0.00024 

u.tlz._oo42 

988 ±31 

0.881^ 
0.850 ±0.035 
1.43 ±0.11 
0.533 ± 0.029 
0.205 ± 0.025 
0.317 ±0.024 
0.275 ± 0.010 


64.8 ± 1.7 
1.080 ±0.057 
0.500 ± 0.035 
0.060 ± 0.013 
167 ± 26 
0.0548 ± 0.0017 
13.65 ±0.46 

88.25 + q"2^ 

11 -0.054 

5652.889 ±0.035 

TQ7+0.054 
KJ - Jy -0.047 

995 ± 34 
0.884 J|« 
0.863+ o34 
1.37 ±0.14 
0.532 ± 0.029 
0.206 ± 0.025 
0.316 ±0.024 
0.275 ± 0.010 


56.8 ± 2.4 
1.115 ±0.047 
0.378 ± 0.027 


0.0379 ±0.0011 
8.35 ±0.16 
84.89 ±0.19 
0.745 ± 0.014 
5584.76697 ± 0.00028 

n 77Q+0.030 
\J.^l J_ o26 

1369 ± 39 
938 +008 ° 

-0.076 

0.976 ± 0.034 
1.0098 ±0.06 
0.424 ± 0.029 
0.277 ±0.017 
0.286 ± 0.020 
0.263 ± 0.020 
0.263 ± 0.020 
0.260 ±0.010 


Radial Velocity RMS 


CORALIE 
HARPS 

RMS/error HARPS 
RMS/error CORALIE 


12.12 
5.76 
1.00 
0.96 


10.7 
4.62 
0.80 
0.85 


13.53 
1.56 



Table 3. Planetary and stellar parameters for WASP-42 and WASP-49, as well as the radial-velocity RMS. For WASP-42 we show 
both the circular and the (preferred) eccentric model. 
(a) p re( iicted 

Assuming A=0 and F=l . 



taken from tables, it might not be accurate for WASP-42. For ex- 
ample, a higher v mac would produce a lower vsin;*, and thus a 
slower rotation and older gyrochronological age. Together with 
the absence of periodic brightness modulations (although we can 
not exclude non-periodic photometric variations of a few milli- 
magnitudes), the evidence points towards an older age of several 
Gyrs for WASP-42. 

The planet's position in the mass-radius diagram (see Figure 
[8]) is at the low-mass end of the Hot Jupiters. We used the tab- 
ulated radii of |Fortney et al.| ( |2007| l in order to check whether 
the planetary radius is within predictions using the values 
for 0. 1 AU = 875 K, and interpolating for a planetary mass of 
0.5 Mj. Assuming the gyrochronologically determined age of 1 
Gyr, we obtain a core mass of slightly above 10 M e , while using 
values for 4.5 Gyrs would suggest a core mass close to M®. 
Compared to other planets of similar mass, WASP-42 b is no 



outlier in terms of density, but with p p = 0.40 + 0.05 p j falls into 
the well-populated region between 0.25 pj and 0.5 pj. We also 
compared the planet's radius to the radius predicted from the fits 
describing the Saturn-mass population presented in |Enoch et al. 
( |2012| l. Using Equation 8 of Eno ch et al.| ( |2012"| ) together with 
our uncertainties on stellar and planet parameters, we obtain a 
predicted radius of 1.31 ± 0.03 Rj, which is within the scatter 
around the model (Figure 14 in Enoch et al. ( 2012| l). 



Adopting an eccentric orbit for WASP-42 b allows for a 
small (0.060 ± 0.013) but non-zero eccentricity, while the argu- 
ment of periastron is found to be a> = 167 + 26 deg. We checked 
any dependence of the derived eccentricity on the jitter values 
assumed by re-running the analysis adding jitters of 3 m s~' and 
15 m s~' quadratically to the CORALIE errors as well as a jit- 
ters of 5 m s~' to the HARPS data. In all cases, the significance 
of the derived eccentricity remained above 3 <x. To illustrate, 



8 



M. Lendl et al.: WASP-42 b and WASP-49 b: two new transiting sub-Jupiters. 



the distribution of the argument of periastron and eccentricity 
found from the analysis including a jitter of 5 m s _1 for HARPS 
and 15 m s _1 for CORALIE is depicted in Figure [9] In order to 
evaluate the robustness of this detection we computed the Bayes 
Factor (e.g. Carlin & Louis 2009) of the two models, resulting 
in a value of B ec = 2440. This is a very strong indication that 
the eccentric model is a better representation of the data. To pro- 
vide a second evaluation we compared the reduced chi-squared 
values xi . of the radial-velocity models, for datasets including 
all data (CORALIE and HARPS), as well as all data but leav- 
ing out the radial-velocity sequence obtained on April 4 2011, 
for and the HARPS or CORALIE data separately. For the lat- 
ter cases we re-ran the MCMC analyses and found eccentricities 
compatible with the one derived from the complete dataset. We 
also used the prescription of |Lucy & Sweeney| ( |1971[ ) (LS test) in 
order to test each of the above for the significance of the circular 
solution. The results are summarized in Table [4] In short, in all 
cases, the non-circular model gives a better fit to the data, and the 
values found for the eccentricity agree from all subsets. The LS 
test indicates a probability of 99.5 % for the eccentric model to 
be required if the HARPS data are considered alone, and a 99 % 
probability if the entire dataset is considered, both above the 5 % 
limit suggested by |Lucy & Sweeney ( 1971| l. The CORALIE data 
alone do not justify the eccentric model. Considering that there 
are only 10 HARPS points outside the consecutive time series, 
more data would be highly desirable in order to secure the detec- 
tion. Based on these arguments, we conclude that there is con- 
siderable evidence for an eccentric orbit of WASP-42 b and we 
present the eccentric as well as the circular solution in Table [3] 

The planet is located on the outer edge of the known Hot- 
Jupiter pileup, at a separation of 4.05 times the Roche limit, well 
above the 2 Rrl cutoff identified by |Ford & Rasio ( |2006| l, be- 
low which planets are thought to be unable to attain that loca- 
tion by the circularization of previously eccentric orbits. This 
means WASP-42 b is compatible with the evolutionary scenario 
described by, for example, Matsu mura et al.| l f2010) : formation at 
farther orbital separations, scattering via planet-planet or Kozai 
interactions to an eccentric orbit, and subsequent tidal circular- 
ization. 



4.2. WASP-49 b 



is a near Saturn-mass (M = 0.38 ± 0.03My) 
a 2. 8 -day orbit around the G6 star 



WASP-49 
planet in 

2MASS 06042146-1657550. In W ASP-49, lack of lithium 
would suggest an age of several Gy ( Sestito & Randich 2005 ) 
and the low vsin/* implies a long rotation period of around 
50 d. The gyrochronological relation of Barnes (2007) gives 
an age of ~ 13+g 5 Gy. Hence, we conclude that WASP-49 is a 
relatively old main-sequence star. We performed an isochrone 
analysis the same way as for WASP-42 with z = 0.011 (Figure 
7(b) 1, finding an age above 10 Gyr, in good agreement with the 



age found from gyrochronology. 

In the period-mass plane WASP-49 b lies clearly below the 
bulk of Hot Jupiters and is nearly a twin of the well known 
HD 149026 b (|Sato et al.||2005]l, a planet known for its high 



density. While the two planets have almost identical mass and 
period, they show very different radii and orbit two very differ- 
ent stars - while WASP-49 is remarkable for having one of the 
lowest metallicities known for planet host stars (-0.23 + 0.07), 
HD 149026 is one of the most metal-rich planet hosts. In contrast 
to HD 149026 b, WASP-49 b has a low density (0.27 ± 0.03 pj). 
WASP-49 b fills a gap between 0.32 and 0.42 My in the mass- 
radius plane (see Figure IS), being the first low-density planet in 




Fig. 8. Mass-Radius diagram of transiting exoplanets. WASP- 
42 b and WASP-49 b are depicted using upward and down- 
ward facing triangles, respectively. The dotted lines represent 
curves of equal density at (from top to bottom) p = 0.25 pj, 
p = 0.5 pj and p — pj. The planet parameters were taken from 
www.exoplanet.eu. 



this mass range. Using the tabulated values from |Fortney et alT 
(2007) for system ages of 4.5 Gyrs, an equilibrium temperature 
of 1300 K and interpolating for a planet mass of 0.38 Mj, we find 
a predicted maximal planet radius of 1 .08 R j (the case of a zero 
core mass). As with WASP-42 b, we employed Equation 8 of 
Enoch et al. (2012) to obtain a predicted radius of 1.26 + 0.02 Rj. 
Again, this value is within the scatter around the fit (Figure 14 in 
Enoch etal.| ( |20T2| ). We conclude that WASP-49 b, while larger 
than predicted from models, is not an outlier in the set of known 
planets. 

There is no evidence for a non-zero eccentricity of WASP- 
49 b, and we determine a 3-cr upper limit of e — 0.09. WASP- 
49 b is separated from its host star by 2.4 times the Roche Limit, 
a rather low value which, together with the planet's old age, 
might indicate that orbital decay has occurred since the arrival 
of the planet at short orbital distances. 

To summarize, WASP-42 b and WASP-49 b are two newly 
discovered close-in transiting planets with masses in the range 
1-2 times that of Saturn. While WASP-49 b is old, inflated and 
orbiting close to a metal-poor G6 star, WASP-42 b is orbiting a 
cooler Kl star at the outer edge of the well-known planet pile-up, 
probably on a slightly eccentric orbit. 
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Dataset e XL.** xLl^O LStest ° 

HARPS + CORALIE 0.060!" !!" 0.82 ±0.18 1.33 ±0.22 0.010 

HARPS + CORALIE (no RM sequence) 0.061!™ 0.87 ± 0.23 1.53 ±0.30 0.024 

CORALIE 0.068!° : °g 1.13 ±0.34 1.33 ±0.34 0.24 

HARPS - 059 ^:o!3 0.90 ±0.26 1.58 ±0.33 0.005 

Table 4. Eccentricities, x 2 values and LS-test results of the eccentric and circular models for WASP-42 b found from all RV data 
and the datasets from the two instruments separately. 

(a) 



Probability of the circular model to be accurate, as defined by Lucy & Sweeney (1971 
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Fig. 9. The distribution of eccentricity and argument of perias- 
tron for WASP-42 b obtained from the global analysis while 
adding a jitter of 5 m s _1 to the HARPS and 15 m s to the 
CORALIE data. 
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